Derivation of Newton's Law of Gravitation Based on a Fluid Mechanical Singularity 

Model of Particles 
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We speculate that the universe may be filled with a kind of fluid which may be called aether or 
tao. Thus, Newton's law of gravitation is derived by methods of hydrodynamics based on a sink 
flow model of particles. 
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I. INTRODUCTION 



Newton's law of gravitation can be written as 



F21 



mim 2 A 
G — 5— r 2 i, 



(1) 



where m\ and m 2 are the masses of two particles, r is 
the distance between the two particles, G is the gravita- 
tional constant, F12 is the force exerted on the particle 
with mass m 2 by the particle with mass mi, f 2 i denotes 
the unit vector directed outward along the line from the 
particle with mass mi to the particle with mass m 2 . 

The main purpose of this paper is to derive Newton's 
law of gravitation by means of fluid mechanics based on 
sink flow model of particles. 

The motive of this paper is to seek a mechanism of 
gravitational phenomena. The reasons why new models 
of gravity are interesting may be summarized as follows. 

Firstly, there exists some astronomical phenomena 
that could not be interpreted by the present theories of 
gravitation, for instance, the Titius-Bode law lj. New 
theories of gravity may view these problems from new 
angles. 

Secondly, whether the gravitational constant G de- 
pends on time and space is still unknown, for instance, 
see HUH 110,011. It is known that the gravitational 
constant G is a constant in Newton's theory of gravita- 
tion and in theory of general relativity. 

Thirdly, the mechanism of the action-at-a-distance 
gravitation remains an unsolved problem in physics for 
more than 300 years [1 0, 0] . Althou gh t heory of gen- 
eral relativity is a field theory of gravity |l2j, the concept 
of field is different from that of continuum mechanics 
fisl Q, Hi! Hfij ] because of the absence of a continuum in 
theory of general relativity. Thus, theory of general rela- 
tivity can only be regarded as a phenomenological theory 
of gravity. 

Fourthly, we do not have a satisfactory quantum theory 
of gravity presently El El El El E3 0ne of the most 
important challenges in theoretical physics is to reconcile 
quantum theory and theory of general relativity [ItI . 



New theories of gravity may open new ways to solve this 
problem. 

Fifthly, one of the puzzles in physics is the problem of 
dark matter and dark energy, for instance, see [H liil 

mm ni mm 113. New theories of gravity may 

provide new methods to attack this problem |24l l25| . 

Finally, we do not have a successful unified field the- 
ory presently. Great progress has been made towards 
an unification of the four fundamental interactions in the 
universe in the 20th century. However, gravitation is still 
not unified successfully. New theories of gravity may shed 
some light on this puzzle. 

To conclude, it seems that new considerations on grav- 
itation is needed. It is worthy keeping an open mind with 
respect to all the theories of gravitaty before the above 
problems been solved. 

Now let us briefly review the long history of me- 
chanical interpretations of gravitational phenomena. 
Many philosophers and scientists, such as |22j, Thales, 
Anaximenes, believed that everything in the universe is 
made of a kind of fundamental substance^]- Descartes 
was the first to bring the concept of aether into sci- 
ence by suggesting that it has mechanical properties!! . 
Descartes interpreted the celestial motions of celestial 
bodies based on the hypothesis that the universe is filled 
by an fluidic vortex aether. Since Newton's law of 
gravitation was published in 1687 [3^, this action-at-a- 
distance theory was criticized by the French CartesianH. 
Newton admitted that his law did not touch on the mech- 
anism of gravitation]^. He tried to obtain a derivation 
of his law based on Descartes' scientific research program. 
At last, he proved that Descartes' vortex aether hypothe- 
sis could not explain celestial motions properly |33j. New- 
ton himself even suggested an explanation of gravity 
based on the action of an aetherial medium pervading 
the spacefill H3- Euler attempted to explain gravitaty 
based on some hypotheses of a fluidic aether 0. 

In a remarkable paper published in 1905, Einstein 
abandoned the concept of aether j3|| . However, Einstein's 
assertion did not cease the explorations of aether, for 
instance, see 0, H3, US EH E3, E3] . Einstein 

changed his view later and introduced his new concept 



1 



of ether |44l l45j. I regret to admit that it is impossible 
for me to mention all the works related to this field in 
history. 

Inspired by the above mentioned thoughts and 
othersIS El El El E3, we show that Newton's law 
of gravitation is derived based on the hypothesis that all 
the particles are made of singularities of a kind of ideal 
fluid. 



II. FORCES ACTING ON SOURCES AND SINKS 
IN IDEAL FLUIDS 



The purpose of this section is to calculate the forces be- 
tween sources and sinks in inviscid incompressible fluids 
which is called ideal fluids usually. 

Suppose the velocity field u of an ideal fluid is irrota- 
tional, then we have 0, El El El El 



u = V</>, 



(2) 



where 4> is the velocity potential, V = i^+j^+k^ is 
the Hamilton operator. 

It is known that the equation of mass conserv ation of 
an ideal fluid becomes Laplace's equation [l^HIT 



El 



\7 2 <f> = 0, 



(3) 



i):r- 



dz 1 



where <f> is velocity potential, V 2 = 
the Laplace operator. 

Using spherical coordinates (r, 9, <p), a general form of 
solution of Laplace's equation (0 can be obtained by 
sepatation of variables asjsfij 



<Kr,6) = T,(Air l + 



1=0 



r l+i 



Pi (cos ( 



(4) 



where Ai and B>i are arbitrary constants, Pi(x) are Leg- 
endre's function of the first kind which is defined as 



(5) 



If there exists a velocity field which is continuous and 
finite at all points of the space, with the exception of 
individual isolated points, then these isolated points are 
called singularities usually. 

Definition 1 Suppose there exists a singularity at point 
Pa = {xq, yo, zq). If the velocity field of the singularity at 
point P = (x, y, z) is 



u(x,y,z,t) 



Anr 2 



r, 



(0) 



where r = y/(x — xq) 2 + (y — yo) 2 + (z — zq) 2 , f denotes 
the unit vector directed outward along the line from the 
singularity to the point P = (x,y,z), then we call this 
singularity a source if Q > or a sink if Q < 0. Q is 
called the strength of the source or sink. 



From (J2J, we see that the velocity potential 4>(r,9) of a 
source or sink is a solution of Laplace's equation V 2 = 0. 
A calculation shows that the strength Q of a source or 
sink evaluates the volume of the fluid leaving or entering 
a control surface per unit time. 

Theorem 2 Suppose (l)there exists an ideal fluid (2)the 
ideal fluid is irrotational and barotropic, (3)the density 
p is homogeneous, that is dp/dx — dp/dy — dp/dz = 
dp/dt — (4)there are no external body forces exerted 
on the fluid, (5)the fluid is unbounded and the velocity of 
the fluid at the infinity is approaching to zero. Suppose a 
source or sink is stationary and is immersed in the ideal 
fluid. Then, there is a force 



Q 



-pQu 



(7) 



exerted on the source by the fluid, where p is the density 
of the fluid, Q is the strength of the source or the sink, 
Uq is the velocity of the fluid at the location of the source 
induced by all means other than the source itself. 

Proof. Only the proof of the case of a source is needed. 
Let us select the coordinates that is attached to the static 
fluid at the infinity. 

We set the origin of the coordinates at the location 
of the source. We surround the source by an arbitrary 
small spherical surface S. The surface S is centered at 
the origin of the coordinates with radius r. The outward 
unit normal to the surface S is denoted by n. Let r(t) 
denotes the mass system of fluid enclosed in the volume 
between the surface S and the source at time t. Let Fq 
denotes the hydrodynamic force exerted on the source 
by the mass system r, then a reaction of this force must 
act on the the fluid enclosed in the mass system r. Let 
Fs denotes the hydrodynamic force exerted on the mass 
system r due to the pressure distribution on the surface 
S, K denotes momentum of the mass system r. 

As an application of Newton's second law of motion to 
the mass system r,we have 



DK 

~Dt 



(8) 



where D/Di represents the material derivative in the la- 
grangian system. The expressions of the momentum K 
and the force Fg are 



K 



pudV, 



(-p)ndS, 



(9) 



where the first integral is volume integral, the second 
integral is surface integral, n denotes the unit vector di- 
rected outward along the line from the origin of the co- 
ordinates to the field point(x, y, z). 

Since the velocity field is irrotational, we have the fol- 
lowing relation 



V0, 



(10) 
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where <fi is the velocity potential. 

According to Ostrogradsky-Gauss theorem, we have 



pudV 



pVcjydV = gj) pcjyndS. (11) 



Note that the mass system r does not include the sin- 
gularity at the origin. According to Reynolds' transport 
theorem, we have 

5 /// p " dv ~ I III padv + #/ u(u ■ n)ds - (12) 

where V is the volume fixed in space which coincide with 
the mass system r(i) at time t, that is V = r(t). 
Then, using lfT2)l , © and (fTTfl . we have 



DK 

TjF 



P-j^ndS + (JJ) pu{u ■ n)dS. (13) 



According to Lagrange-Cauchy integral 52], we have 
d<f> (V0) 2 p 



at 



... ■ - = /(*), 

2 p 



(14) 



where /(f) is an arbitrary function of time f. Since the 
velocity u of the fluid at the infinity is approaching to 
zero, and noticing |@}, </>(f) must be of the following form 



(r,M) = £^#PKcos0), 

(=0 



(15) 



where Bi(t),l > are functions of time £. Thus, we have 
the following estimations at the infinity of the velocity 
field 



(16) 



where y?(a;) = O(-0(ie)), a; — > a stands for lim^^a | ip(x) \ 
/ip(x) = k, (0 < k < +oo). 

Applying (|14H at the infinity and using l|16|) . we have 
u |— » 0, d<f>/dt — > andp = p,^, where Poo is a constant. 
Thus, /(f) = Poo/ P- Therefore, according to l(Ti)l. we 
have 



Using © and l(T7)l. we have 



Using ©, dJl, we have 



p(u • u) 



/j(u • u)n 



d5. 



(17) 



(18) 



Now let us calculate this velocity u in order to obtain 
Fq. Since the velocity field induced by the source Q 
is ©, then according to the superposition principle of 
velocity field of ideal fluids, the velocity on the surface S 



is 



u 



Q 

Aitr 2 



n + u , 



(20) 



where n denotes the unit vector directed outward along 
the line from the origin of the coordinates to the field 
point (a;, y, z). Using (|19fl and l|20() . we have 



= P 



Q 2 



327r 2 r 4 



n 



-(u ■u )n 



4 7rr 2 



u - (uo • n)u 



dS. 



(21) 



Since the radius r can be arbitrarily small, the velocity 
Uo can be treated as a constant in the integral of l|21|l . 
Thus, |j2U turns out to be 



= ~P 



Q 

Aitr 2 



u dS. 



(22) 



Since again Uo can be treated as a constant, 1)22(1 turns 
out to be (JJJl.This completes the proof. □ 

Remark. Lagally 46J, Landweber and Yih[47l lisj ]. 
Faber 49] and Currie 50J obtained the same result of The- 
orem |21 for the special case where the velocity field is 
steady. 

Theorem|21only considers the situation that the sources 
or sinks are at rest. Now let us consider the case that 
the sources or sinks are moving in the fluid. 

Theorem 3 Suppose the presuppositions (1),(2),(3),(4) 
and (5)in Theorem^ are valid and a source or a sink 
is moving in the fluid with a velocity ~v s , then there is a 
force 



Fq = -pQ(u f - V s 



(23) 



is exerted on the source by the fluid, where p is the density 
of the fluid, Q is the strength of the source or the sink, 
Uf is the velocity of the fluid at the location of the source 
induced by all means other than the source itself. 

Proof. The velocity of the fluid relative to the source 
at the location of the source is u/ — v s . Let us select the 
coordinates that is attached to the source and set the 
origin of the coordinates at the location of the source. 
Then can be arrived following the same procedures 
in the proof of Theorem [21 □ 

Applying Theorem |3 to the situation that a source or 
sink is exposed to the velocity field of another source or 
sink, we have 



~p(u • u)n - pu(u ■ n) 



dS. 



(19) 



Corollary 4 Suppose the presuppositions (1),(2),(3),(4) 
and (5) in Theorem^ are valid and a source or a sink 
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with strength Qi is exposed to the velocity field of another 
source or sink with strength Q\, then the force F21 exerted 
on the singularity with strength Q2 by the velocity field 
of the singularity with strength Qi is 



Theorem 7 The eguation of motion of a particle is 



n Ql - 



PQ2V2, 



(24) 



where r 2 i denotes the unit vector directed outward along 
the line from the singularity with strength Qi to the sin- 
gularity with strength Q2, r is the distance between the 
two singularities, v 2 is the velocity of the source with 
strength Q 2 . 



III. DERIVATION OF INVERSE-SQUARE-LAW 
OF GRAVITATION 



Since quantum theory shows that vacuum is not empty 
and has physical effects, e.g., the Casimir effect 0i El 
13 EH]> it is valuable to probe vacuum by introducing the 
following hypotheses 

Hypothesis 5 Suppose the universe is filled by an ideal 
fluid named tao; the ideal fluid fulfill the conditions 
(2), (3), (4), (5) in Theorem^ 

This fluid may be named tao after pi3 ] in order to distin- 
guish with Cartesian aether. Following Einstein, Infold 
and Hoffmann, who introduced the idea that particles 
may be looked as singularities in fields |5fiu57| . and notic- 
ing (|24|) . it is nature to introduce the following 

Hypothesis 6 All the microscopic particles were made 
up of a kind of elementary sinks of tao. These elementary 
sinks were created simultaneously. The initial masses and 
the strengthes of the elementary sinks are the same. 

We may call these elementary sinks as monads. 

Suppose a particle with mass m is composed of N mon- 
ads. Then, according to Hypothesis H3 we have 



m (t) = m (0) + pq t, 

Q = -Nq , m(t) = Nm (t) 



Q 



(25) 

m (t), (26) 



dmo 



dm 



dt =Pq °> 1£ = - PQ ' 



(27) 



where mo(t) is the mass of monad at time t, — qo(qo > 0) 
is the strength of a monad, m(t) is the mass of a particle 
at time t, Q is the strength of the particle, N is the 
number of monads that make up the particle, p is the 
density of the tao, t ^ 0. 

From l|27l) . we see that the mass mo of a monad is 
increasing since q evaluates the volume of the tao fluid 
entering the monad per unit time. From (|27|l . we also 
see that the mass of a monad or a particle is increasing 
linearly. 

Based on Hypothesis |S] and Hypothesis the motion 
of a particle is determined by 



, . dv 



pqo 



dt mo(t) 



m(t)u 



pqo 

m Q (t) 



m(t)v + F, (28) 



where mo(t) is the mass of monad at time t, —qo is the 
strength of a monad, m(t) is the mass of a particle at time 
t, v is the velocity of the particle, u is the velocity of the 
tao at the location of the particle induced by all means 
other than the particle itself, F denotes other forces. 

Proof. Applying Newton's second law and Theorem 
13 to this particle, we have radv/di = —pQ(u — v) + F. 
Noticing p ^l . we get PS || . □ 

128|) shows that there exists an universal damping force 



pqo 
— ? 

to 



(29) 



exerted on each particle. 

Now let us consider a system consists of two particles. 
Based on Hypothesis [SJ applying Theorem [7| to this sys- 
tem, we have 

Corollary 8 Suppose there is a system consists of two 
particles and there are no other forces exerted on the par- 
ticles, then the equations of motion of this system are 



mi 



TO 2 ' 



<ivi 
~dt 
dv 2 
~dt 



pqo r 

m ' 

pqo r 

mo ' 



-TO1V1 



-TO2V2 



pq TO!TO 2 

~. n o Tl2 \6\)) 

pql 171x1712 , ( v 

"j 2 2~ r21 ' \ 6L ) 

47TTOo r z 



where m,i,i = 1,2 is the mass of the particles, Vj,« = 
1,2 is the velocity of the particles, too is the mass of a 
monad, —qo is the strength of a monad, p is the density 
of the tao, i 12 denotes the unit vector directed outward 
along the line from the particle with mass m 2 (t) to the 
particle with mass m\(t), f 2 i denotes the unit vector di- 
rected outward along the line from the particle with mass 
TOi(i) to the particle with mass mi (i) . 

Ignoring the damping forces in (|31|l . we have 

Corollary 9 Suppose (l)vi <C u,-,i = 1,2, where Vj is 
the velocity of the particle with mass mi, Ui is the velocity 
of the tao at the location of the particle with mass mi 
induced by the other particle, (2)there are no other forces 
exerted on the particles, then the force F2i(t) exerted on 
the particle with mass TO2 (t) by the velocity field of tao 
induced by the particle with mass m\ (t) is 



F 21 (t) = -G(t) 



TOi(i)m. 2 (i) 



"*2i, 



(32) 



where G — pq^ / (Airm^t)) , ?2i denotes the unit vector di- 
rected outward along the line from the particle with mass 
TOi (t) to the particle with mass m^i), T is the distance 
between the two particles. 
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Corollary El is coincide with Newton's inverse-square-law 
of gravitation (JIJ except for two differences. The first 
difference is that mi,i = 1,2 are constants in Newton's 
law Q while in Corollary OH they are functions of time t. 
The second difference is that G is a constant in Newton's 
theory while here G is a function of time t. 

Let us now introduce a hypothesis that G and the 
masses of particles are changing so slowly relative to the 
time scaler of human beings that they can be treated as 
constants approximately. Thus, Newton's law of gravita- 
tion may be considered as a result of Corollary El based 
on this hypothesis. 



IV. SUPERPOSITION PRINCIPLE OF 
GRAVITATIONAL FIELD 

The definition of gravitational field G of a particle with 
mass m is G = F /m test , where m tes t is the mass of a test 
point mass, F is the gravitational force exerted on the 
test point mass by the gravitational field of the particle 
with mass m. Based on Theorem and Corollary El wc 
have 



pqo 

T 

m 



-u. 



(33) 



where p is the density of the tao, mo is the mass of a 
monad, go is the strength of a monad, u is the velocity of 
the tao at the location of the test point mass induced by 
the particle with mass m. Thus, the superposition princi- 
ple of gravitational field is nothing but the superposition 
theorem of the velocity field of ideal fluids. 



V. EQUATIONS OF GRAVITATIONAL FIELD 
OF CONTINUOUSLY DISTRIBUTED 
PARTICLES 

The definition of the volume density of continuously 
distributed sink is 



Pa 



urn — — , 

av^o AV 



(34) 



where Ay is a small volume, AQ is the sum of strengthes 
of all the sinks in the volume AV. 

Now let us to derive the continuity equation of the 
tao from the principle of mass conservation. Consider 
an arbitrary volume V bounded by a closed surface S 
and fixed in space. Suppose there are some point sources 
or sinks continuously distributed in the volume V. The 
total mass in volume V is 



M 



pdV, 



(35) 



The rate of increase of the total mass in volume V is 

(36) 

The rate of mass outflow through the surface S is 



<)p(u-n)dS= /// V.(pv)dV, (37) 
is JJJv 

where u is the velocity field of the tao. 

The rate of mass created inside the volume V is 



pp s dV. 



(38) 



Now according to the principle of mass conservation, 
and make use of Eq. fiflft , Eq. and Eq. (JSEJ , we have 



pp s dV - 



V • (pv)dV (39) 



Since the volume V is arbitrary, we have 
dp 



Of 



V • (pv) = pp s 



(40) 



Since the tao is homogeneous, that is §§ = = §f = 
|f = 0, Eq.gOJ becomes 

V • u = p s . (41) 

Thus, from (|33() . we have V • G = pqop 8 /rriQ. Let 
p m denotes the volume mass density of continuously 
distributed particles, then according to Hypothesis El 
Eq.(j2E|)and Eq.(g3}, we have 

m Q p s 



(42) 



Thus, noticing that the velocity field u of tao is irrota- 
tional and G = pq$ / (Airing) , the equations of gravita- 
tional field of continuously distributed particles can be 
summarized as 



V x G = 0, 

V • G = -AnGpr, 



(43) 



VI. 



THE EQUIVALENCE OF INERTIAL MASS 
AND GRAVITATIONAL MASS 



where p is the density of the tao. 



According to Hypothesis El and Corollary El we have 
niinertiai = m grav it a tionah where m inert iai is the inertial 
mass of a particle, m gravitational is the inertial mass of 
the particle. 



VII. TIME DEPENDENCE OF 
GRAVITATIONAL CONSTANT G AND MASS 

According to Hypothesis El and Corollary El we have 
G = pq$ / (Airrn^(t)) , where mo(t) is the mass of monad 
at time t, —qo is the strength of a monad, p is the density 
of the tao. The time dependence of gravitational mass 
can be seen from (I27|l . 
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VIII. CONCLUSION 

There are two features of this theory of gravitation. 
The hrst feature is that the gravitational interactions are 
transmitted by a kind of fiuidic medium. The second fea- 
ture is the time dependence of gravitational constant and 
gravitational mass. Newton's law of gravitation is arrived 
if we introduce a hypothesis that G and the masses of 
particles are changing so slowly that they can be treated 
as constants. Further, the equations of gravitational field 
of continuously distributed particles and the equivalence 
of inertial mass and gravitational mass are derived. As 
a byproduct, it is shown that there exists a universal 
damping force exerted on each particle. 



IX. DISCUSSION 

There exists some interesting theoretical, experimental 
and applied problems in the fields of fluid mechanics, 
astronomy, astrophysics and other related fields involving 
this theory of gravitation. 
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